Spin diffusion in one-dimensional classical Heisenberg model 
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The problem of spin diffusion is studied numerically in one-dimensional classical Heisenberg model 
using a deterministic odd even spin precession dynamics. We demonstrate that spin diffusion in this 
model, like energy diffusion, is normal and one obtains a long time diffusive tail in the decay of 
autocorrelation function (ACF). Some variations of the model with different coupling schemes and 
with anisotropy are also studied and we find normal diffusion in all of them. A systematic finite 
size analysis of the Heisenberg model also suggests diffusive spreading of fiuctuation, contrary to 
previous claims of anomalous diffusion. 



The classical Heisenberg model [l|, Q has been exten- 
sively studied, both analytically and numerically, for sev- 
eral decades and has become a prototypical model for 
magnetic insulators. However, one important question 
that still awaits a conclusive answer is regarding the time 
dependent behavior of the spins, particularly at very high 
temperature. In the hydrodynamic limit, the dominant 
mode of fluctuation spreading in this system is believed 
to obey of the standard diffusion phenomenology. In ab- 
sence of any microscopic theoretical formalism, studies 
were mostly numerical, generally involving calculation of 
time correlation functions. Although thephenomenol- 
ogy of spin diffusion is an old concept d, 0] , its validity 
in classical Heisenberg model has been vigourously de- 
bated in recent times. Although much effort jl3l - l20| has 
been devoted to understand whether spin diffusion in this 
system is normal or anomalous, a convincing conclusion 
is yet to be reached. Settling this question is not only 
conceptually important e.g., in understanding transport 
properties of spin systems, but also has direct implica- 
tions in routinely performed experiments e.g., NMR and 
ESR in magnetic compounds p-Q- In the following, we 
present a brief outline of the diffusion phenomenology 
and review some of the earlier studies in this direction. 

Let us consider a one-dimensional chain containing 
Heisenberg spins {Si} (three dimensional unit vectors) 
where, i = 1,2, . . . , N with periodic boundary conditions, 
Sn+i = Si. The Hamiltonian is given by. 



i.e 



n = - 



N 

E 



Ki Si ■ Si+i, 



(1) 



where Ki is the interaction strength between the spins 
Si and S'i+i; the spin-spin coupling is ferromagnetic for 
Ki > and anti- ferromagnetic if Ki < 0. The micro- 
scopic equation of motion can be written as. 



—rrSi — Si X Bi, 

at 



(2) 
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where Bi = A'i_i5i_i + KiSi+i is the local molecular 
field experienced by the spin at site i. Clearly, Eq. ([2]) 
conserves (i) the total energy E = ^ ■ Ei = — ^ ■ Ki Si ■ 
Si+i, and (ii) the total spin S — Si. 

Since there is no long range order in this system at 
any finite temperature and because of the conservation 
of total spin, the spin fluctuation in the hydrodynamic 
limit is expected to follow a continuity (diffusion) equa- 
tion dtSq{t) = —Dsq^Sq{t), where Sq{t) is the (discrete) 
Fourier transform of Si{t) and Dg is the spin diffusion 
constant. A similar equation holds for the energy den- 
sity (since total energy is also a constant of motion). 
The continuity equation implies that in the hydrody- 
namic limit (small q and large t) the spin-spin corre- 
lation function As{q,t) = {Sq{t) ■ S-q{0)) decays with 
time exponentially, i.e., As{q,t) ^ e~^=^^*. A direct con- 
sequence of this is that the spin autocorrelation function 
(ACF) As{t) = jJJ2i{^^{t) ■ St{0)) at late times decays 
with a power law tail As{t) ^ t^". As predicted by 
the diffusion phenomenology, the exponent a is equal to 
1/2 in one dimension. This is also true for energy ACF 

Ae{t)^^J:^{m)E^{0))-t-^^'. 

Recently, the problem of spin diffusion in this model 
was studied by Miiller [l^l and it was reported that 
a = 0.609 ± 0.005 (ID) in the hydrodynamic limit, thus 
significantly differing from the spin diffusion prediction. 
Following this, Gerling et. al. performed extensive 
numerical studies with larger system sizes and for longer 
times. They strongly opposed the claim made in Ref. 

and demonstrated that the slope of the ACF slowly 
decreases as t is increased. Nevertheless, they suggested 
that the problem is computationally difffcult since the 
non-asymptotic behaviour of the spin ACF is quite pro- 
nounced. In yet another work [l^, it was concluded that 
in numerical simulation it is not possible to observe the 
t~^/^ behavior, even if it exists, due to the fact that the 
numerical scheme introduces computational errors and 
this violated the conservation of total spin S. The error 
propagation affects the decay of the ACF and makes it 
anomalous. It was suggested that the correlation func- 
tion may show a crossover from non-diffusive to diffusive 
behavior and the characteristic crossover time will de- 



2 



pend on the precision of the numerical scheme employed. 

Another subsequent numerical work in this direction 
[l6j however claimed that although energy diffusion is 
normal, spin diffusion has an anomalous behavior. jl7j . 
A coupled-mode theory of spin fluctuation [I^ suggested 
that spin diffusion is anomalous with As(t) ~ t~^^^ 
asymptotically. Authors in Ref. [T^ studied few vari- 
ants of the Heisenberg models (alternate coupling, ran- 
dom coupling etc.) and suggested that spin diffusion is 
probably normal with alternate coupling but is anoma- 
lous with random coupling. A nonlinear dynamics study 
[20j presented numerical results in support of anoma- 
lous diffusion and computed z « 1.67, which implies 
a = z"-^ w 0.6. Thus, most of the previous works re- 
fute the validity of normal spin diffusion in this system. 

In this work, we re-investigate spin diffusion in classi- 
cal Heisenberg model on a ring, using a discrete time odd 
even dynamics (DTOE) [ll|. We compute the temporal 
decay of the energy and spin ACFs, and its cumulative 
average (defined later). We perform extensive simulation 
on large system sizes {N = 5000 — 20000) and for very 
large times {t = 10^). To the best of our knowledge, 
such large scale simulation has never been performed in 
this system to study spin diffusion. We demonstrate con- 
vincingly that spin diffusion is normal and the exponent 
a = 1/2 (within numerical accuracy). We also study 
a few other variants of the usual Heisenberg model and 
show that diffusion process is also normal in those cases. 
For small system sizes, the ACFs saturate after a char- 
acteristic timescale that depends on the system size. We 
have performed systematic finite size analysis for small 
systems which again indicates that spin diffusion is in- 
deed normal. In the following, we describe our numerical 
scheme and present the results. 

a. DTOE dynamics: This method of integrating the 
discretised version of equation of motion has been dis- 
cussed in detail elsewhere (see Ref. [HI). For the sake 
of completeness, we present an outline of the DTOE dy- 
namics here. 

To integrate the equation of motion numerically, one 
would naively consider a finite difference equation of the 
Euler form 



s,. 



At 



SxB 



(3) 



and update all the spins at time t and obtain their values 
at the next time-step t + 1. However, it can be shown 
that using Eq. ^ directly and updating all the spins 
simultaneously, lead to the violation of the conservation 
laws stated above for all_At > 0; the length of the spins is 

A way to naturally preserve 



HI 



also not held constant 
the length of spins \Si\ is by using an alternative spin 
precession update equation (instead of Eq. [3]) 



i,t+l 



Scos4) + {Sx B)siii(l>+{S-B)B{1 - cost 



(4) 

where R = BjlBA and 4 = \BAAt [H. This will, 



however, still violate the conservation of total energy and 
total spin. 

For the conservation of total energy, we use an odd- 
even spin update rule where the dynamics described in 
Eq. @ is numerically implemented by alternate parallel 
updates of the spins on odd and even sublatticcs. Thus, 
at each step, first, only even spins are updated using the 
spin precession dynamics Eq. Q while the odd spins 
are kept unaltered. Next, the spins on the odd sublat- 
tice are similarly updated. It is straight forward to check 
that update of any spin 5*^ affects only the energy of the 
neighbouring bonds e^-i and Ci, but their sum (e^-i -l-e^) 
remains constant. Thus DTOE dynamics conserves en- 
ergy strictly and also naturally maintains the individual 
spin lengths. 

However, the total spin S does not remain conserved 
using DTOE dynamics. This is a general problem with 
any standard integration scheme; the conservations are 
only approximately maintained depending on the accu- 
racy of the scheme. 
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FIG. 1. (Color online) Typical evolution of the magnitude 
of the total spin 5 = |S| and its components Sx,Sy, Sz using 
DTOE dynamics with N = 5000, At = 0.05 and with uniform 
coupling Ki — 1. All the four quantities are scaled by a factor 
of 1/N. 

Wc will show that DTOE dynamics is still a better nu- 
merical scheme compared to other conventional schemes 
(e.g. Euler method, Runge-Kutta method) since the ac- 
curacy here is naturally higher for any given value of 
At (as E and \Si\ remain accurately conserved for any 
arbitrary At). Thus, one can choose a relatively larger 
At without accumulating large numerical errors and thus 
probe the time dependent behavior of the system at very 
late times. Although choosing a larger At essentially con- 
verts the equation of motion to a map, it can be shown 
that independent of the value of At, this dynamics al- 
lows the system to settle to the correct equilibrium state 
[m and one obtains the correct static spin correlations. 
Also, we have verified it thoroughly that the asymptotic 
behavior of the ACFs remains unaltered with smaller At, 
only the computation becomes more time consuming. 

The time evolution of the total spin and its compo- 
nents, from typical run using the DTOE dynamics, is 
shown in Fig. [TJ We find that up to large times (t ^ 10^), 
the total spin conservation is approximately maintained; 
the magnitude S and the components Sx,Sy,Sz do not 
show any trend of an overall growth (or decay) with time. 
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An Eulcr-like scheme for the same values of the param- 
eters N and At will however develop numerical instabil- 
ities and 'blow up' much before t ^ 10^. Thus, using 
DTOE dynamics, one can perform extensive numerical 
simulation of the system and reliably determine the be- 
havior of the system in the hydrodynamic limit. Below, 
we mention the details of our numerical simulation and 
present the results. 

b. Results: First, we simulate the case for which all 
interactions are positive (ferromagnetic) and uniform i.e., 
Ki = K (and set to unity without the loss of general- 
ity). Starting from a random initial spin configuration, 
we evolve the system using DTOE dynamics. We com- 
pute the ACFs Ae (t) and Ag (t) , and study their late time 
power law decay. However, the ACF has a slow conver- 
gence to its asymptotic behavior. Again, in course of 
its decay with time, the ACF either saturates if the sys- 
tem size is small or, if the system size is large enough, 
the ACF continues to decay and its the numerical value 
keeps decreasing. As such, it becomes computationally 
more and more expensive to get well averaged data at 
late times. To circumvent these problems, instead of the 
the ACF, we compute its cumulative average, defined as 



C,(i) = 7 f As{T)d7 

r Jo 



(5) 



where As{t) = ^E»^i('^KO) • S,{t)). This cumulative 
autocorrelation function (CACF) Cs{t) has the same 
asymptotic time dependence as the ACF and therefore 
at late times Cs(t) ^ if spin diffusion is normal. 

Moreover, this has the added advantage that the data for 
Cs (t) is much less noisy than that of Ag (t) and therefore 
its asymptotic time dependence can be computed with 
high accuracy. Likewise, one can define a cumulative av- 
erage Ce{t) for Ae{t). 

The functions Ce{t) and Cs{t) obtained using DTOE 
dynamics are shown in Fig. [2] At late times (~ t > 10^), 
we find that all the curves show a clear convergence to 
(broken lines in the figures). The decay of 

the energy CACF (Fig. [5^) docs not come as a surprise 
since energy diffusion in this model was already known to 
be normal. It has also been recently shown that energy 
transport in this model obeys Fourier's law [l^, [HI for 
any nonzero temperature. However, the decay of the spin 
CACF (Fig. ^jp) with an exponent a = 1/2 for almost 
two decades 10^ — 10^) is quite interesting. In fact, 
this clearly indicates that, contrary to previous claims, 
spin diffusion is normal in the classical Heisenberg spin 
system and spin fluctuation spreads diffusively. 

We have also simulated this model with other coupling 
schemes, namely, with alternate coupling Ki = (—1)*, 
and random coupling {Ki = ±1 assigned randomly) that 
have been studied in Ref . [l^ . The authors of Ref . [l^ 
suggested that spin diffusion appears to be diffusive for 
the alternate coupling case, whereas, for random cou- 
pling, it is probably non-diffusive. However, for both the 
cases we find a clear t^^/^ behavior at late times. This is 




FIG. 2. (Color online) Log-log plot of the autocorrelation 
functions for the uniform case (a) Ce{t) for A^, At = 20000, 0.1 
averaged over 100 independent realisations, (b) C's{t) for 
N,At = (5000,0.05), (10000,0.05) and (20000,0.1), and av- 
eraged over 200, 100, and 25 independent realisations respec- 
tively. The data points are shifted along the j/-axis for better 
visibility. The broken lines have a slope —0.5 




FIG. 3. (Color online) Log-log plot of the autocorrelation 
function Cs (t) (a) for random coupling (circles) and alternate 
coupling (squares) for N = 20000 with At = 0.05 (b) with 
anisotropy 7 — 0.5 (circles), 1.5 (squares) for A'' — 20000 
with At — 0.10. The data is averaged over 100 independent 
realisations and shifted along j/-axis for clarity. The broken 
lines have a slope —0.5. 



shown in Fig. Thus, our numerical results convinc- 
ingly demonstrate that in all the three models namely, 
uniform, alternate and random coupling, the diffusion 
process is not anomalous. 

We have studied spin diffusion in classical Heisenberg 
model with anisotropic coupling in different spin direc- 
tions 



N 



i+li 



(6) 



where 7 is the anisotropy parameter. Spin diffusion with 
anisotropy in the classical limit has been studied in some 
detail recently (2l|-|23t. Our data for Cs{t) with 7 > 1 
and 7 < 1 is shown in Fig. [5]d. The data shows that at 
late times Cs{t) ~ t~^^^ and thus indicates normal spin 
diffusion in the anisotropic model also. 

Using finite size analysis of the ACF, one can have 
an alternative method of estimating the exponent a [20| . 
We work with small system sizes N < 200 Heisenberg 
spins on a ring with uniform coupling Ki — 1. As stated 
earlier, Aij{t) for small system size saturates at some char- 
acteristic relaxation time ts{N). The saturation of ACF 
for different N is shown in Fig. Numerically, we es- 
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FIG. 4. (Color online) (a) As{t) for finite systems N saturate 
at difi^erent characteristic times fs(A*'). (b) Plot of ts{N)/N'^ 
as a function of A'' which roughly becomes constant for A*" > 
100. The inset shows a log-log plot of ts{N) vs. A'^; the best 
fitted straight line for N > 120 gives a slope z = 2.02 ± 0.04. 



timatc ts{N) by the intersection of a power law fit, for 
the linear part of the curve (in logarithmic scale) before it 
saturates, and the straight line y = As(oo), where ^^(oo) 
is the saturation value of Ag{t) at late times. The satura- 
tion As (oo) is numerically computed by averaging Ag (i) 
far away from the saturation point, i.e., for t ^ tg- The 
relaxation timescale ts (N) is related to the system size N 
via the dynamical exponent z as ts ^ . For a diffusive 
process z should be equal to 2, which is indeed the case 
as can be seen from Fig. HJd where we have plotted ts/N"^ 
against iV; the data attains a constancy (approximately) 
for iV > 100. We fit the data for iV > 100 which is shown 
in the inset of Fig. |4}d. The best fitted straight line in 
the log-log plot has a slope z = 2.02 ± 0.04. Hence, we 
have a = z~^ « 0.5, implying normal spin diffusion. The 
value z = 1.67 reported earlier [2^ seems to be due to 
the smaller system sizes studied there. 

To summarize, we have revisited the spin diffusion 



problem in classical Heisenberg spin model in one di- 
mension. We have performed extensive simulation of 
the model using DTOE dynamics that preserves the con- 
servations of the total energy E accurately. Although 
this dynamics is identical to the equation of the motion 
only in the At limit, it equilibrates the system to 
the correct stationary state for any finite At ll|. It 



is thus advantageous here to use a relatively larger At 
and probe the dynamical behavior of the system up to 
late times. By computing the autocorrelation functions 
we show that, similar to energy diffusion, spin diffusion 
in classical one dimensional Heisenberg model is normal 
{a = 1/2), contrary to what has been suggested in some 
of the previous works. We obtain an estimate for the 
dynamical exponent z « 2, which again indicates that 
spin diffusion is normal. The probable reasons as to why 
most of the previous works concluded that diffusion in 
this system is anomalous could be because of (a) small 
scale simulations - both in system size and time, (b) the 
noisy correlation function data, and (c) the accuracy of 
the method used. Our way of simulation and analysis 
take care of most of these issues and produce a clear long 
time diffusive tails for the correlation function. Although 
this dynamics still lacks the strict total spin conservation, 
however unlike conventional integration schemes where 
the errors accumulation is relatively fast, here the to- 
tal spin conservation is approximately preserved allowing 
one to measure the autocorrelation functions up to very 
large time. It remains a challenge to find a suitable dy- 
namics for this model, which will strictly preserves both 
energy and spin conservation. 
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